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ABSTRACT 



We give distance-redshift relations in terms of elliptic integrals for three different 
mass distributions of the Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology. 
These models are dynamically pressure free FLRW on large scales but, due to mass in- 
homogeneities, differ in their optical properties. They are the filled-beam model (stan- 
dard FLRW), the empty-beam model (no mass density exists in the observing beams) 
and the 2/3 filled-beam model. For special fi m - values the elliptic integrals reduce 
to more familiar functions. These new expressions for distance-redshift significantly 
reduce computer evaluation times. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. INTRODUCTION 



As limits on the global cosmological parameters Vt m and A have been refined, Schmidt et al. 
(1998) and Perlmutter et al. (1999), the optical inadequacy of the standard distance-redshift relation 
(D-z) of FLRW has become more apparent. The problem was first recognized long ago by Zel'dovich 
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(1964), Bertotti (1966), and Kantowski (1969) but the lack of relevant data limited its significance. 
Even though the average mass density parameter Q m (along with Hq and A) determines the large 
scale dynamic behavior of the pressure free universe, knowledge of the actual mass inhomogeneity 
is necessary to accurately determine these parameters from most observations. Most observations 
determine Q, m and A by (indirectly) comparing theoretical D-z curves to observed data. However, 
D-z depends on more than the average mass density. It can depend significantly on details of how 
the mass is distributed, i.e., on how inhomogeneous the mass is on the scale of the widths of the 
observing beams. If some significant fraction {pi/ po < 1) of the total mass density is in the form of 
inhomogeneities and is excluded from the lines of sight to the distant objects observed, a modified, 
i.e., a partially filled-beam D-z is required. 

The necessity of taking into account the effect of inhomogeneities on observations is relatively 
easy to understand. Homogeneous matter inside an observing beam of light gravitationally focuses 
the beam much differently than does an equal-mass clump of externally lensing matter. The simplest 
correction for this gravity- light effect requires the introduction of another parameter is, < v < 2, 
which gives the fraction pi/po = v{v + l)/6 of the mass density of the universe removed from the 
observing beams as inhomogeneities. Using v rather than pi/po or some other parameter is dictated 
by the mathematics of special functions. A reduced mass density in an observing beam causes it 
to diverge relative to a standard FLRW beam. For an observed object in such a universe to have 
the standard FLRW angular size it would thus have to be moved to a smaller z; i.e., objects will 
appear less bright than in the standard FLRW universe. A reasonable application of this model 
to SNe la observations takes pj as the galactic contribution to the total mass density po and the 
remaining contribution as a smooth intergalactic medium. Galaxies are easily excluded from SNe la 
foregrounds by selection (intended or not) and if galaxy mass roughly follows light, including their 
mass in pi is appropriate. In the partially filled-beam model where the additional parameter v ^ 
has been introduced, only lensing by mass clumps external to the beam has been neglected. To 
compare individual observations to D-z of this model requires only an occasional lensing correction; 
however, comparison with the standard FLRW D-z (v = 0) model requires a defocusing correction 
for the partially empty-beam of every observation, as well as the occasional lensing correction. If 
only weak and transparent lensing occurs (to the z max being observed) the standard FLRW D-z 
{v = 0) should give the mean D-z curve. Wang (1999) argues that by using flux-averaging the mean 
can be accurately obtained. Kantowski (1998a) and Kantowski (1998b) claims that determining 
cosmological parameters from data compared with the partially filled Hubble curves given here 
is likely to be easier. Beyond selection effects, unknown lensing probabilities can be highly non- 
Gaussian and should make the mean more difficult to observationally determine, i.e., should require 
more data if a given accuracy of the cosmic parameters is to be obtained, Bertotti (1966); Holz & 
Wald (1998); Holz (1998). The down side for partially filled-beam models is that you must select 
against lensing and must determine the additional parameter v. 

In Sec. 2 we outline the procedure required to obtain D-z for partially filled-beam FLRW 
observations and how the result simplifies for the three special cases of v = 0, 1, and 2. In Sec. 3 
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we give the new results for these three special cases. Some concluding remarks are given in Sec. 4 
and in the Appendix we discuss our Fortran implementation of these results. 



2. The Luminosity Distance-redshift Relation 

For models being discussed here (and for most cosmological models), angular or apparent size 
distance is related to luminosity distance by D < (z) = D((z)/(1 + z) 2 . Hence we need to give only 
one or the other, and we have chosen to give luminosity distances. The Dg(z) which accounts for 
a partially depleted mass density in the observing beam but neglects lensing by external masses is 
found by integrating the second order differential equation for the cross sectional area A(z) of an 
observing beam from source (z = z s ) to observer (z = 0), see Kantowski (1998a) for some history 
of this equation: 



(i + z) V 1 + n m z + n A [(i + z)- 2 - 1] x 

+ (3 + ^ (2 ~^W + zf^Mz) = 0. (1) 



^(1 + Z)V1 + «m* + "a[(1 + Z)- 2 - 1] 



The required boundary conditions are: 



VA\ S = 0, 



Stt—— -, (2) 



dz H a (l + z a y 

where (5f2 is the solid angle of the beam at the source and the FLRW value of the Hubble parameter 
at z s is related to the current value Hq at z = by: 

H s = H (l + Z 9 ) y/l + n m Z s + n A [(l + Z s )-2-l]. (3) 

The luminosity distance is then simply related to the area A\ of the beam at the observer by: 

Equation (1) can be put into the form of a Lame equation and its solution has been given in terms 
of Heun functions in Kantowski (1998a). Solutions can also be given in terms of Lame functions 
but neither Heun nor Lame functions are currently available in standard computer libraries. Con- 
sequently, such expressions are not particularly useful for comparison with data, at this time. For 
the special case where A = the Lame functions reduce to associated Legendre functions and these 
expressions are useful. Other special cases also exist as is pointed out in Kantowski (1998a). 

In the next section we give useful expressions for for three special cases where A is arbitrary 
but where the filling parameter v is restricted to values 0,1, and 2. For these three cases we can write 
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Di as an elliptic integral and hence we can give Di in terms of the three fundamental incomplete 
Legendre elliptic integrals F((j>, k), E((f), k), and 11(0, a 2 , k). These functions are universally available 
and these new expressions significantly speed up the evaluation of Di (see the Appendix). Distance- 
redshift for Q,q = 1 can be given in terms of hypergeometric functions, see (21) and (53), or 
associated Legendre functions, see (22) and (54); however, we also give as more complicated 
expressions involving Legendre elliptic integrals, (23) and (55), because these expressions evaluate 
more rapidly using currently available Fortran routines. 

It is not at all clear that the solution of (1) can be written as elliptic integrals for the special 
cases of v = 0, 1 and 2. However, the steps required to arrive at this conclusion can be found 
in Whittaker & Watson (1927) under integral functions for Lame and Matthew equations (see 
especially Sec. 19.53). The authors have carried out the conversion directly for all three cases; 
however, the v = and 2 conversions can be reached by simpler means. The integral for v = 0, 
the standard FLRW filled-beam case, is given in (5) and is well known. The v = 2 (empty-beam) 
integral given in (46) is easy to obtain because the coefficient of y/A vanishes in (1). The first 
integral is trivial and the second is elliptic resulting in (46). For v = 1, the 66% filled-beam model, 
the integral is given in (30); however, no simple way of getting this from (1) seems to exist. 

In Sec. 3. we outline results for all big bang models in the first quadrant of the fl m - plane 
(see Fig. 1), hoping to facilitate their usage. Luminosity distances for the three large open domains 
are given in subsections A, and for the boundaries of these domains in subsections B. 



3. Luminosity Distances as Legendre Elliptic Integrals 
I. v = 0, Completely Filled-Beam Observations (Standard FLRW) 
A. Three Open Big Bang Domains 

Kaufman & Schucking (1971) and Kaufman (1971) gave magnitude-redshift relations for stan- 
dard pressure- free FLRW models as inverse Weierstrass functions and more recently Feige (1992) 
gave comoving distances and light travel times for these models using Legendre elliptic integrals. 
In this section we give simpler and more useful results which are directly comparable with Ed- 
wards (1972) who used Jacobi elliptic functions. The well known and often used integral form for 
luminosity distance in standard FLRW is: 



c 1 + z 

D e (Q m , fi A , v = 0; z) = — , S K 



Vl 01 Jo x/OTzm 



dz 



^(i + z) 2 {i + n mZ )-z{z + 2)n A 



(5) 



which we integrate using Byrd & Friedman (1971) to obtain, 

-g{F(<p z ,k) - F(4> ,k)}}, 



C 1 + 2 

De(n m , fi A , v = 0; z) = — — = , S K 



H y/\l-n \ 



(6) 
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or equivalently using an addition formula for F((f>,k), i.e., F(<f> z ,k) — F(4>o,k) = F(A(/) z ,k) we get: 



D e (0, m ,n A ,iy = 0;z) 



1 + z 



H, 



-g F(A0 z ,k) 



(7) 



The parameter k = (f2o — l)/|^o ~~ 1| is determined by the sign of the 3-curvature and S K [] is one 
of two functions: 

J sinh[ ] : n = —1, 
1 sinf 1 : k = +1. 



S K [ 



Constants g and k depend on the cosmic parameters Q. m & J7a, and F(<p,k) is the incomplete 
Legendre elliptic integral of the first kind. 1 The constants g and k depend on Q m &; Q\ only 
through a combination called b defined by: 



b =-(27/2)- ^ A 



(i- w 



-oo < b < oo, 



(8) 



b < 44> k = -1, 
6 > k = +1. 

The functions 4> z and A0 2 depend on the redshift z and the cosmic parameters £l m k, Q\ (not just 
on the combination b). Domains for the various b values in the fi m - Q\ plane are shown in Fig. 1. 

1. For the two open domains defined by b < and 2 < b, quantities g, k, 4> z , and A0 2 are 
conveniently written in terms of intermediate constants v K , y\ and A defined by: 



v K = 



K(b-l) + y/b(b-2) 



1/3 



v K > 1. 



-1 + k(v k + v K v ) 

2/1 = o : 



A = A(^ m , n A ) = v / yi(3yi + 2) 

Parameters 5 and k are then given by: 



> 1. 



and 



g = g{VL m , VLk) = l/y/A(Q m ,n A ) 
2A + k(1 + 3yi) 



v 2 K + + 1 



1/4 



< 1, 



k 2 = k 2 (f} m ,ft A ) = 



AA 



\ + \9 2 (Vn + V- 1 ) 



< 1. 



(9) 
(10) 

(11) 

(12) 
(13) 
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Functions (j) z , and A<p z are given by: 

<j) z = 0(fi m , Oa; 2) = cos -1 

and 

Acj) z = A<p(Q m ,Q A ;z) = 2 tan -1 



(i + z)n m /\i-n \ + K yi -A 



_(i + z)a m /\i-n \ + K yi + A\ ' 



(14) 



VAy/\i-n \Ji + z[i - (i - fio)^ 1 ^]- 1 



1 + z[l - (1 - ^o)^ 1 ^]- 1 + V(l + ^) 2 (1 + [! m 2) - z(z + 2)ft A 

(15) " 

2. For the domain 0<6<2(=>k = 1) three intermediate parameter yi, yi and 2/3 are 
convenient to use, although none are really necessary. In this domain of b, intermediate parameters 
yi, ?/2 and 2/2 are related to the cosmic parameters VL m Sz JIa through b by: 



.'/. = 3 ( 1 + «> s 



i/2 = - ( -1 - 2 cos 



2/3 = g i "I + COS 



cos _1 (l - 6) 
3 

cos'^l - b) 
3 

cos -1 (l - 6) 



+ \/3 si 



sin 



cos _1 (l - 6) 



, -i<y 2 <-2/3, 
cos _1 (l - 6) 



— \/3 si 



sin 



, < yi < 1/3, 



, -2/3<2/ 3 <0. (16) 



The following expressions are valid only in the lower right part of the £l m — £l\ plane. In the upper 
left domain where b also satisfies < b < 2, expressions can be given, but there a big bang doesn't 
occur. The parameters g and k and functions (fi z and A(f> z needed to evaluate (6) and (7) are: 



g = g(n m ,n A ) = 



Vm - V2 ' 



k 2 = k 2 (n m ,n A ) = 



yi - 2/3 
yi - 2/2 



< l, 



4>(0, m ,Q A ;z) = sin 1 



yi - 2/2 



(l + z)tt m /|l-ft | + 2/i' 



(17) 

(18) 
(19) 



A<f) z 



A0(Q m ,n A ;z) 

Vyi - 2/2 \Zm - fW(i - ^o) - Vv3-( 1 + z)vt m /{\ - n ) 



= 2 tan -1 



y/[ yi - fi m /(i - fi )][y 2 - (i + *)JW(i - n )] + V[yi ^- m] 



, (20) 



where j/i < — > yi means repeat the previous term with y\ and j/2 exchanged. 
B. Boundaries 
1. Qo = ft m + Q A = 1 
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For the spatially flat model ( b — ► ±00) a much simpler expression involving hyper geometric 
functions results: 



c f z 

D e (n m ,n A = 1 - tt m ,v = 0;z) = — (1 + z) I —= 

no Jo V 1 

c 2(1 + z) 



ciz 



+ n m z(3 + 3z + z 2 ) 



12 7 

2F1 [ 7:, 7T; 7; 1 - 



6' 3' 6' 



2-^1 



1 - Q r 



\l + Q m z(3 + 3z + z 2 )}^ 
When f2 m / 1 (21) can be expressed as associated Legendre functions, 

c 2 1 / f T(l/6)(l + z) 



1 2 7 



6' 3' 6' 1 + n m z(3 + 3z + z 2 



(21) 



A(fi m , Qa = 1 - fim, ^ = 0; z) 



ffo 3^(1 -fi ro )]Vi2 



-1/6 



1 



1 



-1/6 / / 1 + f2 m z(3 + 3z + z 2 ) 



- 1/6 VV^W (1+^)^4) -1/6 ^ O m (l + z)3 

If (22) is given in terms of Legendre elliptic integrals the result is more complicated: 

c 1 + z 



.(22) 



Di(tt m , fl A = l- Q m , v = 0; z) 



where 



H (3) 1 /4^(f7-i _ 1)1/6 
c 1 + z 

h~o (3) 1 /4 v /n^(n- 1 - 1)1/6 



1 \/3 

2 + - 



{F(^,k)-F(^,k)} 
F(A^,k) 



0(J2 m ;z) = cos 1 



i + z + (i-V3)(n-i -1)1/3 

l + z + (l + V3)(^ 1 -l)V3 



(23) 

(24) 
(25) 



and 



A(f> z = A(j)(tt m ,tt A ;z) = 2 tan 1 



Z Jv^ m (lMn " l) 1 /3yi + 2; [l + (l/n m _ 1)1/3] 



-1 



1 + z[l + (l/fi m - l)i/3]-i + y/l + ^ m z(3 + 3z + z 2 ) 



(26) 



2. b = 2 



This value of 6 can be identified with "critical" values of the cosmic parameters, Felten & 
Isaacman (1986). We give a result good only for the lower 6 = 2 curve, see (44). These models 
start with a big bang and expand to the the finite Einstein radius at t = 00, see A3(vii-b) in the 



-8- 



appendix of McVittie (1965): 



D e (Q m , n A (n m ), v = 0; z) 



c 1 + z 
Hoy/\l-n \ 



x sin < In 



v/i/3 - fi m /(i - n ) + 1 



y 1/3 - n m /(i - n„) - 1 



Vi/3-(i + «)n m /(i-n )-i 



Vi/3-(i + «)n m /(i-n ) + i 



3. o A = o 

This result is due to Mattig (1958), we include it for completeness: 



D e (n m ,tt A = 0, i/ = 0;z) 



2c 



z + (O m -2) (vT+n^i-i)} 



(27) 



(28) 



4. ft m = 

These are massless big bang models, Q A < 1, discussed by Robertson (1933): 

D e (n m = o, n A , u = o; Z ) = c ^ + ^ {i + z - y/n A + (i + ^) 2 (i - ^a)} 



(29) 



II. v = l, 66% Filled-Beam Observations 
A. Four Open Big Bang Domains 

D e (Q mi n A ,u = l;z) = 



— 2(l+z)Sign[3-O m /(l-O )] 

#0 



[3 - n m (i + z)/(i - n )} [3 - n ro /(i - ^ 



(l-Oo)[36 + ^ A /(l-^o) 3 ] 
'(n^A,*) v / I(1-^o)[36 + ^a/(1-^o) 3 ]| x 

2 [3 - fi m (l + z)/(l - fi )] Vl 1 + *) 2 (1 + ^W) " <z + 2 )^A 



Jo 



where 



S, 



(30) 



(f2 m ,n A ,z) L 



cosh[ 
sinh[ 
sin[ ] 
sinh[ 



b < o& [3 - n m (i + z)/(i - fi )] [3 - n m /(i - n )] < o, 

b < 0k [3 - O m (l + z)/(l - fi )] [3 - - «o)] > 0, 

< b < 486, 
486 < b. 



Only the principal value of the integral (P) is needed and unlike the v = case, this integral takes 
on different forms when evaluated using Legendre elliptic integrals, depending on the value of the 
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parameter b. Parts of the analytic result (31) sometimes diverge even though the total expression 
remains finite. For example when b = 486, i.e., when ^36 + fi^fi^/Cl — ^o) 3 = or equivalently 
yi = 3, a limit must be taken. The resulting Di on this new boundary can be found in II. B. 5 below. 
This new boundary splits the one open domain 2 < b < oo into two parts, see Fig. 2. Consequently, 
the Q m - JIa plane is more complicated for v = 1 than for either v = or v = 2. See Al below for 
additional trouble points that occur. 

1. For the three open domains defined by b < 0, 2 < b < 486, and 486 < b the luminosity 
distance D( takes the form: 



D e (n m ,n A ,u = l;z) 



— 2(l+z)Sign[3-fi m /(l-« )] 
-Ho 



[3 - n m (i + z)/(i - n )} [3 - n ro /(i - ^ 



'(n m ,f} A ,z) 



k V|36 + ^ A /(1-O )3| 
2 v 7 ! [A + k( Vi - 3)] 



(l-Oo)[36 + ^ A /(l-^o) 
[F(^,k)-F(^,,k)] 



+ 4 „ 3) [p n(^, d 2 , k) - p n(0o, « 2 , k)] 1 + /& 



Mvi - 3 ) 



(31) 



where yi,A, k, and Z are defined in (10)-(14) and the additional constant d 2 is: 

^2 = (A + K ( yi -3)) 2 



AAk{ Vi - 3) 



(32) 



H(4>, a 2 ,k) is the incomplete Legendre elliptic integral of the third kind 2 and P H((j), a 2 , k) is the 
principal part of that integral. The function is one of, 



fb 



\ln\{[l+h(z)}[l-h(0)]}/{[l+h(0)][l-h(z)}}\ : b < or 486 < b, 
\ [tan" 1 h(z) - tan- 1 h(0)] : 2<b< 486, 



where h(z) is defined by: 



h(z) = 



VI36 + gpgg - jgyq + z)n m /\i - n Q \ + jjjjj 

(3 - yi)^/[(l + 2)^/11 - fi | - «(1 + yi)/2] 2 - (1 + - 3 yi )/4 



(33) 



Some care has to be taken when using these expressions. Divergences in the function /& necessarily 
occur and cancel divergences in H((f>, a 2 , k). Divergences in also occur which add to divergences in 
n(</>, a 2 , k) and cancel zeros in the multiplicative factor -y/|[3 — f2 m (l + -z)/(l — ^o)] [3 — ^m/(l — ^o) 
of (31). Redshift independent divergences occur when fi m / (1 — $7q) = 3 and when O m (3 — yi)/(l — 



2 Il(0,a 2 ,k) = J" * 1/ — a 2 sin 2 0) \/l — k 2 sin 2 J dep. In arriving at the results for the two-thirds filled beam 
model we discovered that equation 361.54 of Byrd & Friedman (1971) has the two square-root terms interchanged 
for the case a 2 /(a 2 — 1) > k 2 . 
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^o) = 2/i(2yi + 5). These points are plotted in Figure 2. Redshift dependent divergences occur 
at (1 + z) = 3(1 — $7o)/^m and at (1 + z)Q m (3 — yi)/(l — CIq) = J/i (5 + 2yi). These points ap- 
pear in the O m - CI a plane respectively to the left of the Q m /(1 — Qq) = 3 line and between the 
O m (3 — 2/i)/(l — &o) = + 2yi) and b = 486 curves. 

Computer evaluation of (31) can be speeded up by reducing the number of Legendre elliptic 
integrals that must be evaluated. As in (7) we can use the addition formula for F((f>, k), i.e., 
F((f> z ,k) — F((p Q ,k) = F(A(j) z ,k) and an addition formula for U(<fi, a 2 ,k), 3 



1 



U(</> g , a 2 , k) - n(^ , a\ k) = n(A^, a 2 , k) + - 



cr 



log 



2 I/ (a 2 - l)(a 2 - k 2 ) to V 1 - € 



1 + C 



(34) 



where 



sin c/> 2 sin <pQ sin A^ z y / a 2 (a 2 — l)(a 2 — k 2 ) 



(35) 



1 — a 2 sin 2 A<f> z — a 2 sin <f> z sin 0o cos A<p z \J\ — k 2 sin 2 A0 Z 

to cut the number of elliptic functions from four to two. We were not able to simplify this expression 
enough to justify inclusion of a rewritten version of (31). However, it was used in our Fortran 
implementation (see Appendix). 

2. For the open domain defined by < b < 2 the luminosity distance Dg has a somewhat 
simpler form: 



c 2(1 + 3) [3 - n m (i + z)/(i - n Q )] [3 - n m /(i - n )] 

1; Z = — \ / r— „o „ — 77-, „ ^ X 



sm 



y/36 + 0^/(1 - O ) 3 

(3 - yi)Vm - 2/2 
yi - 3 



[36 + 2 t O A /(l - O ) 3 ] 
F(^,k)-F(0 o ,k) 



+ 



n 



yi - 2/2 



n 



yi - 2/2 



(36) 



The constants yi, yi and k, and the function <fi z are as defined in I. A. 2 above [see (16)-(19)]. Just 
as in the previous case, the number of Legendre elliptic functions in (36) can be reduced from four 
to two by using the appropriate addition formulas. For F(<f), k) the formula is always the same, see 
(6) and (7), but because a 2 is negative (34) changes to: 



II((k, a 2 , k) - 11(00, a 2 , k) = n(A0 z , a 2 , k) 



2 V (l-a 2 )(a 2 -k 2 ) 



tan" 1 (0 , 



where 



sin cj) z sin cfio sin Acp z ^a 2 {l — a 2 )(a 2 — k 2 ) 
I — a 2 sin 2 Acp z — a 2 sin 4> z sin 4>q cos Acf) z y/l — k 2 sin 2 Acf) z 



(37) 



(38) 



3 This equation is 116.03 of Byrd & Friedman (1971), corrected for two sign errors. 
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This is 116.02 of Byrd & Friedman (1971) with one sign error corrected. 
B. Boundaries 

1. f^o = O m + = 1 

For these models b — > ±oo and a much simpler expression results: 



D e (n m , n A = l- tt m , v=l;z) 



2(1 + z) 3 / 2 



H >/l — Q r 



sinh 



Vi - n r , 



r 

7o (i + z)a/T 



Hq\/1 — Q r 



-.{1 + zf 



(1 + zVl + JW(3 + 3;z + z 2 ) 
; \ 1/3 



v/1 + O m z(3 + 3z + z 2 ) + VT^ 7 ^ 



i - vr^f^ 



1/3- 



y/1 + O m z(3 + 3z + z 2 ) - VT^a 



(39) 



This result can be given in terms of Legendre elliptic integrals F((f), k) and n(</>, a 2 ,k); however, 
the authors can think of no useful purpose in doing so. 

2. b = 2 

See the description for the v = case in section LB. 2 including (44) for this "critical" value of 



D e (n m , n A (tt m ), u = l;z) 



c (l + z)y/3/2 

H ^/|i -n | (ii) 



V8 



+ 



1-3 



1 - fi. 



1-3 



n m (i + z) 







i-n 







cos (^=\og(h z ) 



+ 1/1-3 



17 r 



1-0 



1 - 3 



{i + z)n r 



l-n 







sin ( -j= log (/i 2 



where h z is defined by: 



1 + ^1/3-^/(1-^) \ /2/3 + (l + z)O m /(l-O ) 



\i + VV3-(i + z )n m /(i-n )y V 

3. A = 

This result was first given by Dyer & Roeder (1973), 



2/3 + O m /(l - O ) 



D e (Q m ,tt A = 0,v = l;z) 
4. n m = 



c 4 
^o30l 



3 j \ 

-fi m - 1 + -O m Z j V 7 ! + 



3 

— r2 m i 



(40) 



(41) 



(42) 
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This result is exactly the same as the v = result (29). If there is no mass in the universe 
then removing 33% of no mass from the beam changes nothing. 

5. b = 486 

This result is equivalent to the b — > 486 limit of (31) but is simpler to use. Because Q,A(Q m ) is 
double valued for b = constant > 2 , two expressions must be given to draw the b = 486 curve, see 
Fig. 2. For the upper part of the curve: 



Q\(0, m ) = 1 — Q m + 3y2/6 Q m cosh 



cosh 



-i 



b/2 (ft" 1 - 1) 



(43) 



where < ft m < 1/(1 — y/2/b). In this expression hyperbolic cosine analytically becomes cosine 
for ft m > 1/(1 + y/2/b). For the lower part of the curve: 



^A(ftm) = 1 - ftm + 3v2/6 Q m cos 



cos 



-i 



b/2 (i - n- 



+ 7T 



(44) 



where 1 < ft m < 1/(1 — y/2/b). The simplified result is: 

c (1 + z) 

D e (Q m , n A (Q m ), v = 1; z) = 
\/33 + 5\ /, 



F 



2V33 



H v/|l-ft |(33)( 3 / 4 ) 



2V33 



ft m (l + z) 
33 + 5A_ jB 



33 + h\ 



2V33 



2V33 



+2 (33) 



(1/4) 



8 + [2 + lW(l- W 



v/3 - ft m /(l - fi ) [3 + \/33 - fi m /(l - fi )] 

^8 + [2 + (i + 3)n m /(i-n )] 2 " 

V3 - (1 + ^)ft m /(l - fto) [3 + V33 - (1 + ^W(l - fio)] 



(45) 



The arguments of the elliptic functions, 4> z and O , can be calculated from (14) using yi = 3 and 
A = \/33- To reduce the number of elliptic functions needed to evaluate (45), addition formulas 
for F(<p,k) and E((j),k) can be used [see (7), (48), and (49)]. The value of A(j) z is given by (15). 



III. v = 2, Empty-Beam Observations 
A. Three Open Big Bang Domains 

c ,o f z dz 



D e (n m ,n A ,u = 2;z) 



(i + zf y/(i + z) 2 {i + n mZ ) - z{z + 2)n A ' 



(46) 



Like the v = 1 case this integral takes on different forms when evaluated in terms of Legendre 
elliptic integrals, depending on the value of the parameter b. 
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1. For the two open domains denned by b < and 2 < b the luminosity distance takes the 



form: 



D e (0, m ,n A ,iy = 2;z 



c (1+z) 



H Q A { 



{A+H!n] (l + z)[(l + :)Q,„/|1 -Qol + .i + K/y,] Q ( „/|.l -Q„| f ,1 + ,■//.! 

(A- Kyi )^/\i-n \ 



2VA 



F(0 2 ,k)-F(^,k) 



+ y/Ay/\l-n \ 



E(</> z ,k) - E(<fo,k) 



(47) 



where yi,A, k, and </> 2 are defined in (10)-(14). 4 Just as with the result for the v = 1 case, i.e., 
(31), the number of Legendre elliptic integrals required to evaluate (47) can be reduced from four 
to two by using addition formulas 116.01 of Byrd & Friedman (1971). The addition formula for 

E(<f>,k) is: 

E((f> z ,k) — £'(0o, k) = E(A(f> z ,k) — k 2 sin <f> z sin <p sin A0 2 . 



(48) 



For this case 



-k sin (f) z sin sin A0 2 



2[2A + K (l + 3yi)] 



[(l + *)fW(l-« )-yi-KA] 



vjg + gWg - »o) - 2/1] [»m/(i - n ) - g 

X [n m /(l - fi ) - 2/1 - «A] [tan(A0 2 /2) + 1/ tan(A0 2 /2)] ' 
where an expression for tan(A^ 2 /2) is given by (15). 

2. For the domain < b < 2 the luminosity distance takes the form: 

D e (n m ,n A ,v = 2;z) = 



(49) 



c (1 + z) 2 

h n A 



"2/3 



yXj + z) 2 {\ + fl m z) - gi + 2)^ A 1 

(i + z)[(i + z)n m /\i - n \ + y3 ] n m /\i - n \ + y 3 



2/2V 7 ! 1 - 



F(cf> z ,k) - F(cf> Q ,k) 



- Vyl - 2/2\/l 1 - ^ol 



E(<f> z ,k) - E(<f>o,-k) 



.,(50) 



V2/1 - 2/2 

where the constants yi, y%, y 3 and k are defined in (16)-(18) but the function <p z is now defined as 



4> z = <j)(n m ,tt A ;z) = sin 



_i (i + z)n m /\i-n \ + y 2 
I (i + z)n m /ji -fi | + 2/3' 



(51) 



For this case the value of A<^> 2 needed to reduce the number of elliptic integrals is the NEGATIVE 
of that given by (20) for the v = case. When the addition formula (48) is used, an additional 



4 E((j>, k) = v/l -k 2 sin 2 # 
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term is contributed to (50) which can be evaluated using, 

^ m (yi-y 3 )|i-^ol ( - 3/2) (yi-y2) ( - 1/2) 



— k sin (f> z sin 4>q sin A(f) z = 



[(i + z)ni/(i - n ) 2 - (2 + z) yi n m /(i - n ) - 2 yi (i + yi )\ 
[ y2 - n m /(i - n )W(i + z)2(i + n mZ ) - z(z + 2)n A [ m - (i + z)n m /(i - n )} } 



[ y3 - n m /(i - n )] 



52) 



^ [y 3 - (1 + - fio)] 
B. Boundaries 

1. r^o = + = 1 

This case is the b — ► ±oo limit of (46) and a simpler expression containing hyper geometric 
functions results: 



De(CL m , n A = l- Q m , v = 2; z) 



(1 + z) v / l + Q m z(3 + 3z + z 2 ) 



+- ft 1 / 3 



[1 + ft m z(3 + 3^ + z 2 )]5/6 

" 2i?1 (^ ; y ;1_ ^ 



5 1 11 

2 r l l7ui "Tr! 



i - a 



6' 3' 6 ' 1 + Q m z(3 + 3z + z 2 ) 

(53) 



When Q m / 1, (53) can be expressed in terms of associated Legendre functions as, 
D e (Q m , tt A = l- n m , v = 2; z) = 



+ i- 



i 



r(5/6) 

(l + z)y/l + n m z(3 + 3z + z 2 ) ' 2V6 



+ 



5/12 



(1 + *) 



1/4 



,-5/6 f /l + O m z(3 + 3z + z 2 ) 



- P 



-5/6 
1 / 6 



yi + ^ m z(3 + 3z + z 2 ) 1/6 On(l + ^) 3 

When Q m / 1, (53) can also be expressed in terms of Legendre elliptic integrals as, 

c (l + zf{ 



(54) 



De(n m ,£l A = 1 - Q m ,u = 2;z) = 



H i-n r 



-(V3 + 1) (O- 1 -!) 



1 1\V3 



^1 + 0„z(3 + 3z + z 2 ) 



(1 + z) [1 + z + (V3 + 1) (Om 1 - 1) 1/3 ] 1 + ( V3 + 1) (^m - 1) 



>V3 



1 i\l/6 



(V3 + l)(3)V* V ^" (n ™ l) 
+(3) 1 / 4 ^(O m 1 -l) 1/6 



F(^,k)-F(^ ,k) 
E^k) - E(<f>o,k) 



(55) 
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where the constant k is given by (24) and the functions 4> z and Acj) z are given respectively by (25) 
and (26). For this case the additional term needed to use the addition formula (48) in (55) is: 

— k 2 sin 4> z sin 0o sin A(f> z 

z 2(3) 3 / 4 (2 + v/3)\/l-fim 



1 + Z + (1 + V3~) (^m " 1) 1/3 1 + (1 + V3) (^m " 1) 



-1 r \V3 







l + ^l + fi m z(3 + 3z + z 2 ) 


} 






+ 2y/l + tt m z(3 + 3z + z 2 )} 



(56) 



2. b = 2 

See the description for the z/ = case in section LB. 2 including (44) for this "critical" value of 



b: 



D e (n m , n A (n m ),u = 2; z) 

c9n m (i + z) 2 ( i R (i + z )n r 



H 2|l-fi | 3 / 2 1 (1 + z) V 3 1-fio 



1 tt r 



3 1 - fin 



+-, — log 



i - n 







1 + Vl/3-(l + *)On/(l-no) / 2 /3 + *V(1 - "o) 



1 + ^1/3-^/(1-^) V 2/3 + (l + z)fi m /(l-fi ) 



(57) 



3. fi A = 

This result was first given by Dyer & Roeder (1972), 

De(Q m ,Q A = 0, v = 2;z) 

\ 2 T sn_ 

In 



c fi m (l + z) 



+ 



fTo 4(1 - fi m ) 3 / 2 
3 



2(1 -n m ) 



1 + VI - / VT + IU - VI - 



and can be rewritten using the identity 



Vi + n m « + 2 ^ Tzrn 



1 - VI - ^m/ VVl + ^m^+ Vl-^r 
VI + !!m^ 



1 - 



(1 + *) 2 



(58) 



sinh 1 



1 i-n m = i ln / yi + fi m z + yr^n; 

fi m (l + z) 2 n Vvi + fi m z- Vl -7 ^ 



When Q m > 1 equation (58) is analytically continued using VI — — ► ±zVfi m — 1, 
simplifies by using, sinh _1 (ix) = isin _1 (x) to give a form containing only real variables. 
fi m = 1 result for all v was given by Dashevskii & Slysh (1966): 



De(tt m = i,fi A = o, z^;z) 



#o> + i) 



+ _ (l+z)(- 



2^2^ 



(59) 

which 
The 

(60) 
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4. n m = o 

This result is exactly the same as the v = and u = 1 result (29). If there is no mass in the 
universe then removing 100% of no mass from the beam removes nothing. 

4. Conclusions 

We have given useful forms for the luminosity distance in three currently relevant cosmologies. 
They are all dynamically FLRW cosmologies in the large but differ in how gravitating matter 
effects optical observations. The models are labeled by an additional parameter v [y = 0, 1, and 
2) beyond the familiar Ho,Q m , and A. The v = model is standard FLRW where all matter is 
homogeneous and transparent on the scale of the observing beam widths. This model is called 
the 'filled-beam' model. The v = 2 model assumes the opposite; all matter is inhomogeneous 
and excluded from the observing beams. This extreme case is called the 'empty-beam' model. The 
v = 1 model assumes that 1/3 of the mass density of the universe is excluded from observing beams 
and hence it is the 'two-thirds filled-beam' model. These three cases were singled out because their 
distance-redshift relations can be given in terms of incomplete elliptic integrals; functions which 
are universally available in computer libraries and very efficiently evaluated. 5 For the v = 1 and 2 
cases, somewhat simpler expressions than what we have given exist, but only for complex arguments 
of the elliptic integrals. We chose to give expressions whose arguments are real and which can be 
rapidly evaluated. Results are available for all < v < 2 but only in terms of the less familiar and 
unavailable Heun functions, Kantowski (1998a). We have extended the flat space, = 1, results 
given here to arbitrary filling parameter v. These new results will be available shortly. Related 
results have been independantly found by Damianski et al. (2000). A calculation similar to the 
v = 1 case given here is that of the age of the Universe as a function of redshift and can be found 
in Thomas & Kantowski (2000). 

R. Kantowski wishes to thank VP for Research, E. Smith, for funds to support J.K. Kao's 
visit to OU during the summer of 1998 when the first elliptic integral results were obtained. R. 
C. Thomas thanks P. Helbig for discussions of his code, see Kayser et al. (1997), and E. Baron for 
benchmarking discussions. 

A. Appendix 

One expected practical use of the results given in this paper is to speedup distance evaluations 
for the v = 0,1,2 partially filled beam FLRW models. We have implemented and made publicly avail- 



5 The results appearing in Section 3 have been coded and are posted at http:/ /www. nhn.ou.edu~thomas/z2dl. html. 
This code is discussed in the Appendix and compared to the numerical integration times of Kayser et al. (1997). 



-17- 



able a Fortran 90 version of this work called Z2DL (see http://www.nhn.ou.edu~thomas/z2dl.html 
for Z2DL with documentation and extensive CPU-time benchmark results). Z2DL uses Carlson 
elliptic integrals (see Press et al. (1994) and references therein) and results in a fast distance cal- 
culator. We have benchmarked Z2DL by comparing it with the commonly used and fast numerical 
integration routine ANGSIZ (see Kayser et al. (1997)). For a given (Cl m ,U\), the total CPU-time 
required to convert 5 x 10 5 redshifts (equally spaced between z=0 and z=5) to luminosity distance 
using Z2DL and ANGSIZ separately were recorded. By calculating the ratio of ANGSIZ CPU-time 
to Z2DL CPU-time on a grid of points in (Cl m ,Cl\) we have generated three speedup surfaces, one 
for each value of v = 0,1,2 (see Fig. 3 for the u = surface). The results for all three comparisons 
are given as contour plots at the web site. Using an IBM AIX 375 MHz Power III approximately 7 
hours was required to generate each (O m , S1a) grid of 30 x 30 points (minus models without a big 
bang). 

For the purpose of a clearer presentation, we omitted speedup points along the Q m = and 
SI A = lines. Along these boundaries speedup factors are greater than 100. The large open domains 
of the Q m -rt\ plane, i.e., subsection 'A' cases, constitute the majority of models in the grid and 
also those with the least impressive speedup. However, even for these cases, the improvement is 
substantial: typically 17-20 for v = (standard filled beam FLRW), 6-8 for v = 1 (66% filled beam 
FLRW), and 11-13 for u = 2 (empty beam FLRW). 

To gauge the level of agreement between distances computed by ANGSIZ and Z2DL, a finer 
grid of (£l m ,Qk) with 3000 x 3000 points (between and 3 in both directions, also excluding 
models without a big bang) was used. For each (Q, m ,Q\), both routines were used to compute 
luminosity distance for z=l. Most often the results agree to within one part in 10 6 . Cases where 
disagreements greater than one part in 10 3 occur are near the upper b=2 line (see Fig. 1). We 
found that ANGSIZ was giving less accurate distances near this boundary of non-big bang models 
as ANGSIZ documentation explains. 
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Fig. 1. — The Q m -Q^_ plane showing various b domains that require different expressions for 
distance-redshift for all three cases: v = 0,1,2 i.e., filled-beam, 66% filled-beam, and empty- 
beam. 

Fig. 2. — Additional domains in the Q m -fi\ plane for v = 1, i.e., for 66% filled-beam observations, 
where complications due to divergent terms occur in the analytic results. For £l m - ft, a values on 
the dashed and dot-dashed lines, define respectively by £l\ = 1 — fi m 4/3 and f2 m (3 — yi)/(l — Oo) = 
y±(2yi +5), expression (31) must be evaluated by taking a numerical limit. For points to the left of 
the straight dashed line and points between the dot-dashed and b = 486 curves, a single value of z 
exits for which (31) also diverges. These z values are defined respectively by (z + 1) = 3(1 — £lo)/Q, m 
and (l+z)fi m (3— yi)/(l — fio) = J/i(5+2yi). For Q m , and z satisfying either equation a limiting 
process must be used to evaluate Dg via (31), see the Appendix. For points on the divergent b = 486 
curve an analytic limit was obtained in (45). 

Fig. 3. — Contour plot of the £l m -£l\ plane showing speedup factors for Z2DL over ANGSIZ when 
v = (standard filled beam FLRW cosmology). Speedup factors for the other two cases considered 
in this paper, v = 1, 2 i.e., the 66% filled-beam and empty-beam can be found at the web site. 



